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ABSTRACT 

Though observationally rare, damped Lya absorption systems dominate the mass 
density of neutral gas in the Universe. Eleven high redshift damped Lya systems 
covering 2.8 < z < 4.4 were discovered in 26 QSOs from the APM z>4 QSO Survey, 
extending these absorption system surveys to the highest redshifts currently possible. 
Combining our new data set with previous surveys we find that the cosmological mass 
density in neutral gas, Q g , does not rise as steeply prior to z^2 as indicated by previous 
studies. There is evidence in the observed O g for a flattening at z^2 and a possible 
turnover at z~3. When combined with the decline at z > 3.5 in number density per 
unit redshift of damped systems with column densities log iVni> 21 atoms cm~ 2 , these 
results point to an epoch at z > 3 prior to which the highest column density damped 
systems are still forming. We find that over the redshift range 2 < z < 4 the total 
mass in neutral gas is marginally comparable with the total visible mass in stars in 
present day galaxies. However, if one considers the total mass visible in stellar disks 
alone, i.e. excluding galactic bulges, the two values are comparable. We are observing 
a mass of neutral gas comparable to the mass of visible disk stars. Lanzetta, Wolfe & 
Turnshek found that Q(z « 3.5) was twice Q(z w 2), implying a much larger amount 
of star formation must have taken place between z=3.5 and z=2 than is indicated 
by metallicity studies. This created a 'cosmic G-dwarf problem'. The more gradual 
evolution of Q g we find alleviates this. These results have profound implications for 
theories of galaxy formation. 

Key words: cosmology — galaxies: evolution — galaxies: formation — quasars: absorp- 
tion lines 



While the baryonic content of spiral galaxies that are ob- 
served in the present epoch is concentrated in stars, in the 
past this must have been in the form of gas. The principal 
gaseous component in spiral galaxies is HI which has led 
to surveys for absorption systems detected by the damped 
lines they produce (Wolfe et al. 1986 [WTSC]; Lanzetta et al. 
1991 [LWTLMH] ; Lanzetta, Wolfe, & Turnshek 1995 [LWT], 
Wolfe et al. 1995). Damped Lya absorption systems have 
neutral hydrogen column densities of Nui> 2 x 10 20 atoms 
cm -2 and they dominate the baryonic mass contributed 
by HI. We extend the earlier work on damped Lya sys- 
tems to higher redshifts using twenty-six QSOs from the 
APM Damped Lya Survey (Storrie-Lombardi et al. 1996 
[SMIH], Storrie-Lombardi, Irwin & McMahon 1996 [SIM]), 
with eleven candidate or confirmed damped Lya absorption 



1 INTRODUCTION 



systems covering the redshift range 2.8 < z < 4.4 (8 with z 
> 3.5). These data more than triple the redshift path sur- 
veyed at z > 3 and allow the first systematic study up to z 



= 4.7. 



2 EVOLUTION OF Q G ~ BARYONS IN 
NEUTRAL GAS 



The mean cosmological mass density contributed by Lya 
absorbers can be estimated as 




(1) 



giving the current mass density in units of the current crit- 
ical density (LWTLMH). /i is the mean molecular weight of 
the gas which is taken to be 1.3 (75% H and 25% He by 
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mass), ttih is the mass of the hydrogen atom, p cr it is the 
current critical mass density, iV m i n is the low end of the HI 
column density range being investigated, and f(N, z) is the 
column density distribution function. Unfortunately f(N, z) 
is not a simple function and its evolution with redshift is 
difficult to accurately quantify (LWT, SIM). The integral in 
equation [j] can be estimated using 



Nf(N, z)dN 



AX 



(2) 



where AX is the absorption distance interval. The absorp- 
tion distance X is used to remove the redshift dependence 
in the sample and put everything on a comoving coordi- 
nate scale. If the population of absorbers is nonevolving (i.e. 
the product of their space density multiplied by their cross- 
section does not change with redshift) they have a constant 
number density per unit absorption distance. In a standard 
Friedmann universe X is defined as 



X{z) 



§[(l + z) 3/2 -l] ifqo = 0.5; 
Ul + zf-l] ifq = 0. 



(3) 



(Bahcall & Peebles 1969; cf. Tytler 1987). The errors in 
Q g are also difficult to estimate without knowing f(N,z). 
LWTLMH used the standard error in the distribution of Ahi 
which yields zero error if all the column densities in a bin 
are the same. We have estimated the fractional variance in 
f2 9 by comparing the observed distribution of f(N,z) with 
the equivalent Poisson sampling process. This gives 



( a, 



V P 



N, 



(4) 



and l/y^p fractional errors if all the column densities in- 
cluded in a bin are equal. To address uncertainties in f(N, z) 
we also calculated a maximum likelihood estimate of the er- 
rors in the HI column density. We used the power law with 
an exponential turnover form of the column density distri- 
bution function, i.e. the gamma-distribution from SIM 



f(N, z) = (J./N.)(N/N. 



-P e -N/N, 



(•») 



with log N, = 21.63±0.35, /3 = 1.48±0.30, and /„ = 1.77 x 
1(T 2 . Unlike a pure power law this form has a finite integral 
mass. The maximum-likelihood estimates of the errors agree 
well with the fractional variance. 

LWT found that f2 3 inferred from studies of damped 
systems rises with increasing redshift for 0.008 < z < 3.5. For 
the range 3.0 < z < 3.5 which included 4 damped systems 
and was the the highest redshift bin in the study, Q(z « 3.5) 
was twice Q(z ~ 2). This implied a much larger amount of 
star formation must have taken place between z=3.5 and 
z=2 than is indicated by metallicity studies. Specifically, 
Pettini et al. (1994) measured a low mean metallicity in 
damped Lya systems at z«2.2 inferring that they are ob- 
served prior to the bulk of star formation in the disk. The 
LWT result also implied that the bulk of stars in nearby 
galaxies should be metal poor whereas only a small fraction 
of disk stars in the solar neighbourhood are metal poor. This 
result presented a 'cosmic G-dwarf problem' similar to the 
'G-dwarf problem' described in Schmidt (1963) that comes 
about if bright and faint stars formed at the same rate and 
there was no accretion on to the disk of the Galaxy. They 
concluded that the characteristic epoch of metal produc- 



tion in galaxies occurred after the characteristic epoch of 
star formation and that damped Lya absorbers might trace 
disk as well as spheroid evolution. Fall & Pei (1993; and Pei 
& Fall 1995) have argued that obscuration caused by dust 
in damped Lya systems could lead to significant underesti- 
mates of the neutral gas fraction at all redshifts, particularly 
in the range 1 < z < 2, also possibly explaining these results. 

To further investigate these issues, the high redshift 
confirmed and candidate damped Lya systems discovered 
in the APM Damped Lya Survey (SMIH, SIM) have been 
combined with lower redshift samples (WTSC, LWTLMH, 
LWT) to study the evolution of £l g over the range 0.008 < 
z < 4.7. The combined data set includes 44 absorbers in 
366 quasars. The results are tabulated in table 1 and shown 
in figure 1(a) for qo=0 and 1(b) qo=0.5 (Ho=50) follow- 
ing the format presented in Wolfe et al. (1995). The solid 
bins include the entire data set and the dashed bins exclude 
the new high redshift APM dataQ. The region Q s tar is the 
±1(T range for the mass density in stars in nearby galaxies 
(Gnedin & Ostriker 1992). The point at z=0 is the value 
inferred from 21 cm emission from local galaxies (Fall & Pei 
1993; Rao & Briggs 1993). The most striking result is Q g 
does not rise as steeply prior to z~2 as indicated by previ- 
ous studies. There is now evidence for a flattening in Q g at 
z~2 and a possible turnover at z~3. This result, combined 
with the decline at z>3.5 in number density per unit red- 
shift of damped systems with column densities log JVhi> 21 
atoms cm~ 2 (SIM) points to an epoch at z > 3 prior to 
which the largest damped systems are still forming. The 
decrease in number density at high redshift of the highest 
column density absorbers is in marked contrast to the more 
numerous lower column density systems, e.g. Lyman-limit 
systems (A^hi~ 10 18 atoms cm -2 ) (Storrie-Lombardi et al. 
1994; Stengler-Larrea et al. 1995) and Lya fo rest absorbers 



(JVhi~ 10 ld -10 15 cm" 2 ) (Williger et al. 1994) 



The inclusion of the APM survey data for z > 3 reduces 
significantly the value previously found for £l g in the bin 3 
< z < 3.5. Only one damped system is added to the exist- 
ing four in this redshift range, but the absorption distance 
is doubled, which comes in to the calculation of fl g in the 
denominator in equation ^ (see table 2). The additional red- 
shift path added by the APM survey is shown graphically 
by the sensitivity function in figure 6 of SMIH. We find that 
over the redshift range 2 < z < 4 the total mass in neu- 
tral gas (fig) is marginally comparable with the total visible 
mass in stars in present day galaxies (fl 3 t ar ) for qo=0.5. 
However, if one considers the total mass visible in stellar 
disks alone, i.e. excluding galactic bulges, the two values 
are comparable. Using the result from Schechter & Dressier 
(1987) that galactic disks and bulges contribute equally to 
the mass density of the Universe, we are observing a mass of 
neutral gas comparable to the mass of visible disk stars, i.e. 
Sl 9 ~ Qdisk-stars- We note that the uncertainty in the total 
mass in visible stars in the local Universe is comparable with 
our estimates of the mass in neutral gas at z > 2. Given this, 
and the fact that we do not know if damped systems are the 
precursors to galactic disks, bulges, or both, these results 
are difficult to interpret. If we make a plausible correction 



* Excluding the new APM data effectively yields the data set 
analyzed in LWT. 
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for obscuration by dust as advocated by Pei & Fall (1995) 
the Q, g points shown in figure 1(b) would migrate to the po- 
sitions of the open circles shown in figure 2. More work is 
needed to determine the severity of dust obscuration in op- 
tically selected QSO surveys. An estimated 20% correction 
for the neutral gas not in damped Lya systems is shown by 
the open squares. 
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3 IMPLICATIONS FOR GALAXY 
FORMATION THEORIES 

The shape of the f2 s curve has been used by numerous 
authors to constra in theories of galaxy formation and cos- 



Kauffmann & Char- 



mologica l models (Klypin et al. 1995 
lot 1994 [Ma fc Bertschingcr 1994 M° & Miraldc-Escudc 
1994 ). The y nave found that cold+hot dark matter (CHDM) 



models are incompatible with the previous results from the 
damped Lya systems at z~3 as they predict too few galac- 
tic halos. These models need to be reevaluated now that 
the value of Q g (z=3) indicated from ou r new observation s 



has changed significantly. For example, (Klypin et al. 199 



show a peak in Q g at z=3 for various CHDM models. Though 
the error bars on Q g are still far too large to accurately dis- 
criminate between details of cosmological models, the overall 
shape has implications for structure formation. 



4 SUMMARY 

Combining the new data from the APM high redshift survey 
extends studies of the cosmological mass density in neutral 
gas, f2 9 , to z=4.7. We find evidence for a flattening in Q g at 
z~2 and a possible turnover at z~3. Though the turnover is 
not formally significant, when combined with the decline at 
z > 3.5 in number density of damped systems with column 
densities log Nm> 21 atoms cm -2 , these results point to an 
epoch prior to which the largest damped systems are still 
forming. 

Previous studies indicated that Q(a w 3.5) was twice 
Q(z « 2), implying a much larger amount of star formation 
must have taken place between z=3.5 and z=2 than is in- 
dicated by metallicity studies. The more gradual evolution 
we find in Q g alleviates this problem. The results are con- 
sistent with observations of damped Lya systems at z«2.2 
that show a low mean metallicity, suggesting that they are 
observed prior to the bulk of star formation in the disk. We 



Bahcall J.N., Peebles P.J.E. 1969, ApJ, 156, L7 
Fall, S.M. & Pei, Y.C. 1993, ApJ, 402, 479 
Gnedin, N.Y. & Ostriker, J. P. 1992, ApJ, 400,1 
Kauffmann, G., Chariot, S. 1994, ApJ, 430, L97 
Klypin, A., Borgani, S., Holtzman, J., Primack, J. 1995 ApJ, 444, 
1 

Lanzetta, K.M., Wolfe, A.M., Turnshek, D.A., Lu, L., McMahon, 
R.G. & Hazard, C. 1991, ApJS, 77, 1 [LWTLMH] 

Lanzetta, K.M., Wolfe, A.M. & Turnshek, D.A., 1995, ApJ, 440, 
435 [LWT] 

Ma, CP., Bertschinger, E. 1994, ApJ, 434, L25 
Mo, H.J., Miralde-Escude, J. 1994, ApJ, 430, L25 
Pei, Y.C. & Fall, S.M. 1995, ApJ, 454, 69 

Pettini, M., Smith, L.J., Hunstead, RW. & King, D.L. 1994, ApJ, 
426, 79 

Rao, S. & Briggs, F. 1993, ApJ, 419, 515 
Schechter, P.L., & Dressier, A. 1987, AJ, 94, 563 
Schmidt, M. 1963, ApJ, 137, 758 

Stengler-Larrea E.A., et al., Boksenberg A., Steidel C.C., Sargent 
W.L.W., Bahcall J.N., Bergeron J., Hartig G.F., Januzzi B.T., 
Kirhakos S., Savage B.D., Schneider D.P., Turnshek D.A., 
Weymann R.J., 1995, ApJ, 444, 64 

Storrie-Lombardi, L.J., McMahon, R.G., Irwin, M.J., & Hazard, 
C. 1994, ApJ , 427, L13 

Storrie-Lombardi, L.J., McMahon, R.G., Irwin, M.J. & Hazard, 
C. 1996, ApJ , 468, 128 [SMIH] 

Storrie-Lombardi, L.J., Irwin, M.J. & McMahon, R.G. 1996, MN- 
RAS, in press [SIM] 

Tytler D. 1987, ApJ, 321, 49 

Wolfe, A.M., Lanzetta, K.M., Foltz, C.B. & Chaffee, F.H. 1995, 
ApJ, 454, 698 

Wolfe, A.M., Turnshek, D.A., Smith, H.E. & Cohen R.D. 1986, 
ApJS, 61, 249 [WTSC] 

Williger, G.M., Baldwin, J.A., Carswell, R.F., Cooke, A.J., Haz- 
ard, C, Irwin, M.J., McMahon, R.G., & Storrie-Lombardi, 
L.J. 1994, ApJ , 428, 574 



for bias in optically selected quasar samples due to possi- 
ble obscuration by dust in foreground absorbers. Theories 
of galaxy formation and constraints on cosmological models 
utilising £l g should be reevaluated in light of the new obser- 
vational results presented here. The error bars are still very 
large at high redshift and to differentiate between a peak 
versus a flattening in the Q g curve, larger samples of bright 
z > 3.5 quasars are needed to discover damped Lya systems 
with z > 3. 
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Table 1. Data for Figures 



Bin 


DLA 




#of 




qo = o 


qo 


= 0.5 


Rcdshift 





Az 


DLA QSO 


AX 


n g 


AX 


n 3 

[xh 5() 10 3 ] 


Range 






in bin 




[xft 50 10 3 ] 




.008-1.5 


0.64 


47.8 


4 186 


73.1 


0.56±0.32 


58.5 


0.70±0.40 


1.5-2.0 


1.89 


27.9 


4 126 


79.5 


1.21±0.71 


47.1 


2.05±1.19 


2.0-3.0 


2.40 


120.2 


22 176 


415.9 


1.50±0.49 


223.4 


2.79±0.91 


3.0-3.5 


3.17 


24.3 


5 82 


102.0 


1.48±0.72 


49.8 


3.04±1.48 


3.5-4.7 


4.01 


19.2 


9 32 


93.8 


0.85±0.34 


42.5 


1.87±0.75 


Dashed bins, excluding high rcdshift data. 








2.0-3.0 


2.38 


114.6 


21 154 


394.8 


1.56±0.52 


212.6 


2.90±0.96 


3.0-3.5 


3.19 


11.8 


4 56 


48.9 


2.79±1.47 


24.0 


5.68±3.00 



Table 2. Rcdshift Path and Absorption Distance 



Data Set Az AX AX 

(90 = 0) (qo = 0.5) 



3.0 < z < 3.5 



APM Damped Lya Survey 


12.5 


53.1 


25.8 


WTSC + LWTLMH + LWT 


11.8 


48.9 


24.0 


Combined 


24.3 


102.0 


49.8 


z > 3 








APM Damped Lyct Survey 


30.5 


141.2 


65.5 


WTSC + LWTLMH + LWT 


13.0 


54.6 


26.7 


Combined 


43.5 


195.8 


92.2 


< z < 4.7 








APM Damped Lya Survey 


36.1 


162.3 


76.5 


WTSC + LWTLMH + LWT 


203.4 


602.1 


344.8 


Combined 


239.5 


764.4 


421.3 



Figure 1. The mean cosmological mass density in neutral gas, 
Hg, contributed by damped Lyct absorbers for 0.008<z<4.7 for 
(a) qo=0 and (b) qo=0.5. The solid bins include the combined 
data set and the dashed bins exclude the new APM high redshift 
data. The region U s tar is the ±l<x range for the mass density in 
stars in nearby galaxies (Gnedin & Ostriker 1992). The point at 
z=0 is the value inferred from 21 cm emission from local galaxies 
(Fall & Pei 1993; Rao & Briggs 1993). These results arc tabulated 
in table 1. 



Figure 2. The mean cosmological mass density in neutral gas, 
Q g , contributed by damped Lya absorbers for 0.008<z<4.7. The 
solid bins are the combined data set shown in figure 1(b) for 
qo = 0.5. The circles show the observed data points corrected 
for possible dust obscuration using values determined from the 
closed-box/outflow models shown in figure 4(b) of Pei & Fall 
(1995). The squares add an estimated 20% correction for neutral 
gas not in damped Lyo absorbers. The region Q s tar is the ±1ct 
range for the mass density in stars in nearby galaxies (Gnedin & 
Ostriker 1992). The point at z=0 is the value inferred from 21 
cm emission from local galaxies (Fall & Pei 1993; Rao & Briggs 
1993). 
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